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Abstract 

We try to design a simple model exhibiting self-organized criticality, which 
is amenable to a rigorous mathematical analysis. To this end, we modify 
the generalized Ising Curie-Weiss model by implementing an automatic 
control of the inverse temperature. For a class of symmetric distributions 
whose density satisfies some integrability conditions, we prove that the 
sum Sn of the random variables behaves as in the typical critical gene- 
ralized Ising Curie- Weiss model. The fluctuations are of order n^^* and 
the limiting law is C exp( — Aa;'*) dx where C and A are suitable positive 
constants. 
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1 Introduction 

In their famous article [2] , Per Bak, Chao Tang and Kurt Wiesenfeld showed that 

certain complex systems are naturally attracted by critical points, without any 
external intervention. These systems exhibit the phenomenon of self-organized 
criticality. 

Sclf-organizcd criticality can be observed empirically or simulated on a computer 
in various models. However the mathematical analysis of these models turns out 
to be extremely difficult. Even models whose definition is seemingly simple, such 
as those describing the dynamics of a sandpile, are poorly understood. Other 
challenging models are the models for forest fires [14], which are built with the 
help of the classical percolation process. Some simple models of evolutions also 
lead to critical behaviours [6]. 

Our goal here is to design a model exhibiting self-organized criticality, which is as 
simple as possible, and which is amenable to a rigorous mathematical analysis. 
The most widely studied model in statistical mechanics, which exhibits a phase 
transition and presents critical states, is the Ising model. Its mean field version 
is called the Ising Curie- Weiss model (see for instance [9]). It has been extended 
to real-valued spins by Richard S. Ellis and Charles M. Newman [10], in the 
so called generalized Ising Curie- Weiss model. This model is our starting point 
and we will modify it in order to build a system of interacting random variables, 
which exhibits a phenomenon of self-organized criticality. 

Let us first recall the definition and some results on the generalized Ising Curie- 
Weiss model. Let p be a symmetric probability measure on M with positive 
variance and such that 



The generalized Ising Curie- Weiss model associated to p and the inverse tem- 
perature /3 > is defined through an infinite triangular array of real-valued 
random variables (X^)i<fc<n such that, for all n > 1, has the 

distribution 



where Zn{/3) is a normalization. For any n > 1, we set Sn = -|- • • • -|- X". 

When p = {S^i + Si)/2. we recover the classical Ising Curie- Weiss model. 

We denote by L the Log-Laplacc of p (sec section 6). Richard S. Ellis and Thco- 
dor Eisele have shown in [8] that, if L^^\t) < for any t > 0, then there exists a 
map m which is null on ]0, l/u^], real analytic and positive on ]1/(t^, -|-oo[, and 
such that 



The point 1/cr^ is a critical value and the function m cannot be extended ana- 
lytically around l/cr^. The main theorem of [10] states that, if /3 < l/cr^, then, 
under lJ,n,p,0, 






if /3 < l/a2 
if /3>1/ct2 
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Moreover, for (i = there exists k G N\{0, 1} and A > such that, under 

where Ck^\ is a normalization. This is a consequence of [10] and some properties 
of m explained in [8] implying that the function s i — s- L{s^/P) — has a 
unique maximum at whenever j3 < (see [12] for the details). 

Wc will transform the previous probability distribution in order to obtain a 
model which presents a phenomenon of self-organized criticality, i.e., a model 
which evolves towards the; critical state /3 = l/o"^ of the previous model. More 
precisely, the critical generalized Ising Curie- Weiss model is the model where 
{X^, . . . , X") has the distribution 



{xi-\ \-Xnf 

2na2 



i=l 



We search an automatic control of the inverse temperatiirc /5, which would be 
a function of the random variables in the model, so that, when n goes to +oo, 
/3 converges towards the critical value of the model. We start with the following 
observation : if is a sequence of independent random variables with 

identical distribution p, then, by the law of large numbers, 

— — — > a.s 

Jl n— >cx) 

Thus we are tempted to « replace /3 by n {x\ H h a;^)~^ » in the distribution 

^exp(^- — - — jn^^(-^) 

Hence the model we consider in this paper is given by the distribution 

■ 1 (xi H h a;„)^ 



JJdp(a;i) 

2—1 



These considerations suggest that this model should evolve spontaneously to- 
wards a critical state. We will prove rigorously that our model indeed exhibits 
a phenomenon of self-organized criticality. 

Our main result states that, if p has an even density satisfying some integrability 
conditions, then, asymptotically, the sum Sn of the random variables behaves 
as in the typical critical generalized Ising Curie- Weiss model : if p,^ denotes the 
fourth moment of p, then 

In section 2 we define properly our model. We state our main results and the 
strategy for proving them in section 3. Next we split the proofs in the remaining 
sections (4-10). 
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2 The model 

Let p be a probability measure on M, which is not the Dirac mass at 0. We 
consider an infinite triangular array of real- valued random variables (X*)i<i;<„ 
such that for all n > 1, {X^, . . . , X") has the distribution /I„,p, where 

rf/i„,p(a;i,...,a;„) = ^exp I -— j— 1 1{^2+...+^2 >o} [[dp{xi) 

with 

Z„ = / exp Qi^±-ll±^) >o} n '^^(^') 

JR" 3;;^ i- • • • -|- a;„ / 

We define 5„ = + • • • + X;^ and T„ = (Xi)^ + • • • + {Xl^f. 

The indicator function in the density of the distribution helps to avoid 
any problem of definition if p({0}) is positive, since, if p({0}) > 0, the event 

{x\ H h a;^ = 0} may occur with positive probability. We notice that, unlike 

the generalized model, our model is defined for any probability measure. Indeed 
X I — > is a convex function, therefore 
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i=l / \i=l / i=l 



Thus for any n > 1, 1 < Z„ < e"/^ < +00. 

If we choose p = + Si)/2, we obtain the classical Ising Curie- Weiss model 
at the critical value. 



3 Convergence theorems 

We state here our main results. 

By the classical law of large numbers, if p is centered and has variance cr^, then, 
under p®", (Sn/n,Tn/n) converges in probability towards (0,cr^). The next 
theorem shows that, under the law ^n,p> given certain conditions, {Sn/n, T„/n) 
also converges in probability to (0,cr^) 

Theorem 1. Let p be a symmetric probability measure on M. with positive va- 
riance and such that the function 

A : {u, v)^\n [ 6"^+"^' dp{z) 

JTS. 

is finite in an open neighbourhood of (0, 0). We suppose that one of the following 
conditions holds : 

(a) p has a density 

(b) p is the sum of a finite number of Dirac masses 

(c) There exists c > such that p(]0,c[) = 

(d) p({0}) <l/^e 

Then, under pn^p, {S„/n,Tn/n) converges in probability towards {0,a^). 
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By the classical central limit theorem, imder p*^", Sn/\/n converges in distri- 
bution to a normal distribution with mean zero and variance . The following 
theorem, shows that, given certain conditions, under /U„,p, Sn/rfi^^ converges 
towards a specific distribution. 

Theorem 2. Let p he a probability measure on R with a density f satisfying : 

(a) f is even 

(b) There exists vq > such that 



J 



e"°^ f{z) dz < +00 
(c) There exists p 2] such that 

I F (x + y) r {y)\x\'-P dx dy < +00 
Let be the variance of p and let pi be the fourth moment of p. We have 



The convergence can equivalently be rewritten as 



n3/4 

We prove this convergence in section 10. 

The following corollary is a version of theorem 2 with an hypothesis which is 
weaker but easier to check. 

Corollary 3. Let p be a probability measure on M with an even and bounded 
density f such that 

3vo>0 [ e""^^ dp{z) < +00 
Let be the variance of p and let p4 be the fourth moment of p. Then 



Proof. Wc check that the hypothesis of the corollary imply the condition (c) of 
theorem 2. We have 

/ F/\x + y)f'/'iy)\x\-'/'dxdy 

•''([-1,1]^)= 

< ( sup I/I ) / dxdy+ [ f'/\x + y)f^/\y) dx dy 

V[-2,2] / J[-i,i]2 |a;|V^ J{hiA]^r 
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The second inequality is obtained by applying Fubini's theorem. The first term 
is finite and the second too because, by the Cauchy-Schwarz inequality, 

/ \f{x)\'^^dx<\\f\\^ [ f{x)y^dx=\\f\\^ [ e^°=^'/^f{x)'/^e-^°^'/^dx 
Jm Jr Jr 

< ll/lloo e-'^'^'fix) dx'^ ( 6-"°-' dx^ < +00 



Thus, with p = 3/2 e]l,2], the function {x,y) 1 — y fP{x + y)P{y)\x\^~'' is 
integrable. □ 

For instance, if p has a bounded support and a density which is even and conti- 
nuous on it, then the hypothesis of the theorem are fulfilled. 

We end this section by explaining the strategy for proving these results. 

We denote by i^p the law of {Z, Z^) where Z is a random variable with distribu- 
tion p. By proposition 4, a possible approach to obtain a limit law for {Sn,Tn), 
correctly renormalized, under ]ln,p, is to compute the density of i^*" for n large 
enough, when p has a specific density. We will use this approach in the section 5 
in the case of Gaussian distributions. 

A more robust approach to obtain a limit law for (S'„,r„), corrcictly renormali- 
zed, under Jln,p, is to use the theory of large deviations. We denote by p the 
law of {Sn/n,T„/n) under p'^" and by 9n^p the law of (Sn/n,Tn/n) under ]ln,p- 
Proposition 4, presented in the next section, states that, if A is a subset of M^, 
then 

^'n.^C^) =Mn,p((^,^) = -^^exp(^^) l{y>o}cfi7„_p(a;,j/) 



exp y— j l{j/>o} dun,p{x,y) 

I exp l{3/>0} dA>n,p{x, y) 

By convexity of 1 1 — \ t"^, we have < nTn for any n > 1. We define 

A = {(a;,?/)eM2:x2<t/} and A* = A\{(0,0)} 
Hence we have yn,p (^'^) = 0. Therefore 



For n > 1, under p®". 



n ^ n / n 



i=l 

where (X„,X^)„>i is a sequence of independent and identically distributed 
random variables with distribution Vp. Cramer's theorem implies that (t'n,p)n>i 
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satisfies a weak large deviations principle with speed n, governed by the rate 
function 

I:{x,y)i — > sup {xu + yv ~ A{u,v)) 

(u,i;)eM2 

where for any {u,v) £ R^, 

A(w, v) = ln [ 6""+"* dup{s, t) = ln [ 6"^+"^' dp{z) 

We note that A{u, v) can be equal to +00. If we suppose that the function A 
is finite in an open neighbourhood of (0,0), then / is a good rate function and 
{'^n,p)n>i satisfies a large deviations principle with speed n, governed by I (see 
the section 19 of [5] for these results). 

Here is a classical heuristic : as n goes to +00, the distribution dn,p concentrates 
exponentially fast on the minima of the function 

G = / - F - inf (7 - F) 

where F denotes the map {x, y) \ — > x^ / {2y). Thus, if G has a unique minimum 
at (xo,yo) G A*, then, under Jln.p, {Sn/n,Tn/n) converges in probability to 
{xo,yo)- Moreover, for n large enough, can roughly be approximated by the 
distribution G„ exp(— n/(a;, y)) dx dy where G„ is a renormalization constant. 
Thus, for each bounded continuous function h and a, /3 > 0, 



%.(K^)) 



h{{x — xo)n") exp (— nG(x, y)) dx dy 



A* 



/ exp{—nG{x,y)) dxdy 
J A* 

h{x) exp (— nG (xn"" + xq, yn~^ + yo)) dx dy 



/ exp (— nG (xn "+xo,yn ^ + yo)) dxdy 

J A* 



We use then Laplace method. The key point is the study of the function G in 

the neighbourhood of its minimum (xo,yo)- We have to find four values a G N, 
& S N, A > and B > such that, uniformly on a neighbourhood of {xo,yo), 

-nG(xn-'' +xo,yn-^ + yo) -Ax"" - Byf' 

After computing the law of (S'„/n, T„/n) in section 4 and giving some general 
results on the Cramer transform in section 6, we study the minima of / — F 
in section 7. Next we state a variant of Varadhan's lemma in section 8, which 
helps us to prove theorems 1 and 2 respectively in sections 9 and 10. 
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4 Computation of the law of {Sn/n,Tn/n) 

In this section we compute the laws of {Sn,Tn) and {Sn/n,Tn/n) under 'jln,p- 

Proposition 4. We denote by Vp the law of {Z, Z'^) where Z is a random 
variable with the distribution p. Under 'jln,p, the law of {Sn,Tn) is 

We denote by T'n,p the law of {Sn/n,Tn/n) under p®". Under Jln,p, the law of 
{Sn/n,Tn/n) is 

1 / nx^ \ 

z^^^'^K^) ^^y>°^ dvn,p{^, y) 



Proof. Let / : M? — > R be a bounded measurable function. We have 
1 



We define ^ 

h : {x, y) e I — )■ f{x, y) exp ^v>o} 



We have then 



If " 

^iir.MiSn,Tn)) = — / h{Xi + --- + Xn,xl + --- + xl)\\dp{Xi) 

If " 

= ir h{{xi,xl) + --- + {xn,xl))\\dp{xi) 

If " 

= — / h{zi-\ h2n)]T(ifc'p(2;i) 

= ^ [ h{z)du;-{z) 

Hence the announced law of {Sn,T„), under Jln,p- Moreover, we have for any 
(a;i,...,a;„) eM", 

{xi-\ \-Xnf {{xi-\ ha;„)/n)2 



xl + ---+xl " {xl + ---+xl)/n 
Hence 

^'^"^ (H^' ^) ) = ^ L (ir) ^^^>°> 

This ends the proof of the proposition. □ 
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5 The Gaussian case 

In this section, we prove theorem 2 when p is the Normal law AA(0, cr^) with mean 

and variance fj^. Wc use the method of residue to compute the characteristic 
function of and a Fourier inversion formula to get its density. We finish the 
proof with Laplace method. 

For simplicity, we assume that a"^ — 1. We just write v*^ for f*" and we denote 
by its characteristic function. For {u,v) e M^, 



We need some preliminary results. 

The Gamma distribution with shape k > and scale 9 > 0, denoted by r(A:, 6), 
is the probability distribution with density function 



T{k) 6"= 



-IcOO 



with respect to the Lebesgue measure on R, where F denotes the gamma function 
defined by 

r+oo 

yz>0 T{z)= / x'-'^e-'^dx 
Jo 

For A; > and ^ > 0, the characteristic function of the distribution r{k, 6) is 

u e M ^ (1 - eiu)-^ 

The complex logarithm function (or the principle value of complex logarithm), 
denoted by Log, is defined on O = C\] — oo, 0] by 



yz = x + iy€n Log(2;) = - ln(x +y ) +2i arctan 



2 \x + -\/a;2 + y2 

If a e C and z G Cl, then the a-exponentiation of z is defined by 

z" = exp(aLog(2;)) 
We can now prove the following key lemma : 
Lemma 5. LetteR and C e C such that fHe(C) > 0. Then 



I 



Proof. Let t e M and C = a + i6 e C such that fHe(C) > 0. We define 

K{t, = f e""-^^'/2 

We factorize : 

1, 2 W 1./ tb taV t"^ 
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Thus 

The change of variables y = x — tb/\C,\ gives us 

e'"'^^K{t,Q= I e-^^y-'"'/\^^^"'^dy = - lim / e'^^'/^d^ 

fl-)-+oo 

where the last integral is the contour integral of the entire function z i — > e~^^^/^, 
along the segment 71 in the complex plane with end points R + ita/\(^\ and 

-R + if.a/\C\. 

Let 7 be the rectangle in the complex plane joining successively the points 
R + ita/\(^\, —R + ita/\Q\, —R and R. We apply the residue theorem : 

J y 

since z 1 — > exp{—(z^/2) has no pole. We denote 71, 72, 73 and 74 the successive 
edges of the rectangle 7. 

t 

-R + ita/\(,\ -1 R + ita/\C\ 



72 



74 



73 



R 



J73 J-R Jq 

We make the change of variables y = x^ on ]0, +00 [ : 



+00 



<-"/^dx 



2/ e-i^'/'dx= [ e-<y/'^=[ 
Jo Jo yV Jo 



^-iby/2^-ay/2 ^ 



dy 

5X-V2 



since we recognize, up to a normalization factor, the characteristic function of 
the Gamma distribution with shape 1/2 and scale 2/a. Moreover we have 



a\t\ ( aR"^ Ratbx a ( atx\ , , 

a\t\ f aR"^ Ra\tb\ a f at' 
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Likewise 

/ e-^^'/'dz 
Letting R go to +00, we conclude that 

Since r(l/2) = ^/^T, we obtain the identity stated in the lemma. □ 

Proposition 6. If p = J^{0, 1) then the characteristic function of the dis- 
tribution V™ is 

g ^ exp (-| + Log(l - 2it;))) 

Proof. Let (u, v) G M^. Setting ( = 1 - 2iv G {z G C : d\c{z) > }, we have 
^r^{u,v) = {Mu,v)r = ( [ e»-+* 



(2 

Lemma 5 implies that 

*»(■*••" = (^-■' (-5(1^) (' - ^"'""f 

and the proposition is proved. □ 

Once we know the characteristic function of the law z/*", a Fourier inversion 
formula gives us its density. Wc first have to check that is integrable with 
respect to the Lcbcsgue measure on M^. 

Let {u,v) G Since (1 - 2iv)-^ = (1 + 2to)/(1 + 4i;2), we have 
It follows that 

/ \^^{u,v)\dudv= I exp(- T"' ^\ {\^^v'')-^IUudv 



ijR V 

(l + 4z;2)-(»-2)/4rf^ 



27r 

n 



where we used Fubini's theorem in the third integral. The function 

„^ (l + 4^2)-(„-2)/4 

is continuous on M and integrable in the neighbourhood of +00 and —00 if and 
only if n > 4. 
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Proposition 7. If p = Af{0, 1) and n > 5 then f*" has the density 

{x,y)^^^^[^^W^V\^^y^ exp(-|) (^2/- 1^2<„^ 
with respect to the Lebesgue measure on . 

Proof. We have seen that, if n > 5, then is integrable on M^. The Fourier 
inversion formula implies that i/*" has the density 

/„ : {x,y) ^ / e-"^-'y^<^n{u,v)dudv 

with respect to the Lebesgue measure on M?. Let {x,y) G M^. By Fubini's 
theorem, 

fn{x, y) = -7—--:^ / 7" ^. / exp -ixu - — — ] du] dv 



(27r)2 7„ (1 - 2i«)"/2 V7j, ^\ 2{l-2iv) 
K\—x,- — ^^-— \ dv 



(27r)2 7g (1 - 2w)"/2 V 'l-2w. 
where K is defined by 

Va > V(t, 6) e K{t, a + ib)= [ e'*^-(«+*)^'/2 rf^ 

Lemma 5 implies that for any u € M, 



27r / x2(l-2w)\ /l + 4t;2y/^ 



exp 



/ \ l + 2iv) 



n \ 2n + 
27 / a;2(l-2w)\ 

^"H- 2n J (1-2-)^/^ 



n 



Thus 



'2'Kn 

Therefore \j2'Kn cxp(a:;2/2n)/„(a;, y) is the inverse Fourier transform of the dis- 
tribution r((n — l)/2, 2) taken at the point y — jn. Hence 

„2n(»-3)/2 



y 

n 



( y 
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Finally 

/„(^,y) = (V2^r(^)2("-i)/^) (y-^) exp(-|) 

The proposition is proved. □ 

This previous result and proposition 4 imply that, for n > 5, under the 
law of (S'n.T^) on IR2 is 

Cn^exp(^— - -J (^y- — J 1^2^„ydxdy 
where C„ = Z„ \/2^r((n - l)/2). 

Let a,/? g]0, 1], n > 5 and / a bounded measurable function. The change of 
variables {x, y) i — > (n"a;, n^y) yields 

v[^'^)) = ^iJ^^'y^'^^[-^y T 

2a-l 2\("-3)/2. 



X 



[n'^y-n x ) In2a^2^ni3+iydxdy 



Factorizing by n*^" '^^Z^, we notice that all the terms in the integral are functions 
of x^/n^~^" and y/n^^^. We obtain the following proposition. 

Proposition 8. Let a, /3 g]0, 1]. If p = Af{0, 1) and n > 5 then, under lin,p, 
the distribution of {Sn/n",Tn/n^) is 



„a+/3„(n-3)/2 / / ^2 y \\ / 



Cn 



-^1 



^ ^[z^^ZT^ dxdy 



where X is the indicator function of the set 

£^ = { (x, y) e : y > a; > } 

and ^ and ip are the functions defined on by 

ip:{x,y)< — ^ ^ (-- + 2/ - - x)^ 



(p : {x,y) I — > (y-x) 



-3/2 



We give next some properties of the map ijj in order to determine which values 
of a and (3 to choose. 

Lemma 9. The map tp has a unique minimum at (0, 1) and tp{0, 1) = 1/2. The 
map tp is on and it satisfies : 
■k In the neighbourhood of (0, 1), 

^(a;, y)-\ = \{^^ + {y- 1)') + o{\\x, y - If) 
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★ There exists d > such that for all {x, y) G , 

\x\<5,\y-l\<5 =^ ^(a;,y)-i>i(x2 + (y-l)2) 

★ inf { : |x| > 5 or \y-l\>S}>l/2 

The map if is bounded by 1, it converges to 1 when {x,y) goes to (0, 1) and 
f e-2'^(^''^'V(a;^y)lx=<^;(^a;c^^/ < +oo 

Proof. The map ip is on and, for fixed y > 0, 

dilj, , 1 f 1 1 \ „ 
d^^'^^y^ = 2[-y + i^x)^' 

Equality holds if and only if x = 0. Thus x i — > tp{x,y) is increasing on ]Q,y[ 
and ■0(0)2/) — iv ~ lii(2/))/2- Hence for any (x, y) G , 

^x, y)>^{y- Hy)) > ^ = V(0, 1) 

with equality if and only if (a;, y) = (0, 1). Therefore (0, 1) is the unique minimum 
of ip. In the neighbourhood of (0,0), 

^{x, l + h) = ^{-xil -h + oih^)) + l + h- {h-x-^ih~ xf + o{{h - xf) 

Thus in the neighbourhood of (0, 1), 

i^{x,y)-l = \{x' + {y-ir) + oi\\x,y-ir) 

It follows that there exists S > such that for {x,y) G D^, if \x\ < S and 
|y — 1| < 5, then, 

^ix,y)-^>^ix' + {y-lf) 
Moreover, ii \y — 1\> S and x G [0,y[, then 

tp{x,y) > imin{l-<5-ln(l-(5), l + S-\n{l + 6)} > ^ 
We suppose that 5 <1, otherwise we reduce S. U x > S and y > x, then 

2ij{x,y) >--+y- Hy -S)> inf ( -- +y- ln(y - ^) ) > 1 

y y>s \ y / 

since 6 0. Therefore 

U3f{tp{x,y) :\x\>S or |y - 1| > <5 } > 1/2 
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Finally the map is bounded by 1, it converges to 1 when {x,y) goes to (0, 1) 
and 

f e-^'l'^''^'y^ip{x'^,y)%2^ydxdy <e( f e"^' dx) ( ^ dy) < +oo 

Jr^ ^Jr ^ \Jo Vy ) 

where we used the change of variables (x, y) i — \ {x,y — x"^). □ 
Thus, for fixed {x,y), when n goes to +oo, 

Hence we take a = 3/4 and /3 = 1. We prove now the following theorem : 
Theorem 10. For p = J\f{Q,l), under fin, p, 

Sn c e-^'l^dx . c , 

and — — \ 1 



j7,3/4 „^+oo J^e~y*/'^ dy n n^+oo 



Proof. Let n € N and let / : M.^ — > M be a continuous bounded function. By 

proposition 8, we have 



Let 5 be as in lemma 9. We denote 

„2 



An= I I f{x,y)eii-p{-ni){^,y]] ip{^,y]l^^^2dxdy 

The change of variables (a;, y) i — > (a;, yj \fn + 1) gives 

V^e"/2^„ = / / / ^ + l) exp (-n^ ^ + l)) 



Lemma 9 and the continuity of / imply that 

\s/n s/n ) 2 n^-+cx) 4 4 

Moreover the function inside the integral is dominated by 



(a;,y)^||/||ooexp(-l(x^+y2)) 
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which is independent of n and integrable with respect to the Lebesgue measure 
on R^. By Lebesgue's dominated convergence theorem, we have 

^/^e'^/'^An — ^ / f{x,l)e-''^/^e-y^/^dxdy = V4^[f{x,l)e-''^/^dx 

We define 

Bs = {{x,y)eD^ : \x\ <S,\y-l\<S} 

and 

Bn= [ f{x,y)exp(-ml}(^,y)) (p (^,y) 1^^^2 dxdy 

Let £ = inf { il){x, y) : (x, y) € B|}, 

\Bn\ < e-("-2)^||/||oo j^^exp(-2V>(^,2/)) <^(-^,y) t^^^^.dxdy 
The change of variables {x,y) i — > (xn^/^,j/) yields 

^A^e"/2|B„| < e2^||/|Ue-"(^-V2)n3/4 / e-^^(-'^y^v{x\y)l^.^ydxdy 

Lemma 9 guarantees that e > 1/2 and that the above integrable is finite. The- 
refore s/ne^f^Bn goes to as n goes to +oo. Finally 

j^^ f{x,y)exp (^-ni; (^^,y^^ if (-^'2/) '^V^y>x^ dxdy = An + 

J /(a;,l)e-^*/^da; + o(l)+o(l)) 



e-"/2 



+00 y/n 
^ / fix,l)e--"/Ux 



If / = 1, we have 



+00 

Hence 



/■ / e^^''/"' dx 

-- / /(a;, 2/) -7^ 



"^1(2/) 



This ends the proof of the theorem. □ 
A straightforward change of variables implies that, if p = AA(0, cr^), then, under 

Sn c e-^'/^^"dx , C 2 

and — — > a 



n3/4 n-)-+oo e-y^/*'^'^ dy n n^+x, 
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We also get 

^ A 

Since ^4 = 3cr^, we have 

12a^ ~ 40-4 

We have thus proved theorems 1 and 2 for a Gaussian distribution j\A(0,o-^). 

In the general case we cannot compute explicitly the distribution In the 
following sections, we deal with the more robust method we suggested in the 
heuristics of section 3. 



6 General results on the Cramer transform 

This section, which may be omitted on a first reading, presents some general 
results on the Cramer transform of a probability distribution in M**. Let u he a 
probability measure on R'^, d> 1. The Log-Laplace L of is defined in R"^ by 

VA e L{X) = In / exp(A, z) dv{z) 

where ( , ) denotes the inner product in R''. The Log-Laplace L is convex on R'' 
and takes its values in ] — 00, +00]. We denote by Dl the set where L is finite. 
In particular, if v has a bounded support, then Dl = R"^. 

o 

The^set Dl is convex and contains since L(0) = 0. If Dl 7^ then L is C°° 
on Dl and 

VA eDL,ae N'^ ^^IP^(X) = f ^^"1 . . . exp(A, z) dv{z) 

We refer to [7] and [9] for the proofs of these results. 
We define the Cramer transform of v by 

JixeW^^ sup ( (A, x) - L{X) ) 

It is the Fenchel-Legendre transform of L. Wc write 

Dj = {x eR'^ : J{x) < +00 } 

Proposition 11. (a) J is a non-negative convex and lower semi- continuous 
function. 

(b) If L is finite in a neighbourhood of the origin then the level sets of J, 

{x E R'' : J{x) < a}, a E R, are compact. 

(c) IfteDL and u = VL{t) then ,]{u) = {t,u) - L{t). 

(d) If u has a finite first moment then J{m) = where 




Moreover, ifO€ Dl, then J has a unique minimum at m. 
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Proof. We refer to [7] for the proof of the points (a), (b) and (c). We prove the 
point (d) (see for instance chapter V. of [15]). Let A G R'^. Jensen's inequahty 
implies that 



e<^'^> di^{x) > exp [ (A, x) dv{x) = e^^-™) 



Therefore L{X) > (A, m) and thus J{m) < 0. Since J is a non-negative function, 
it follows that J{m) ~ hence m is a minimum of J. We show that it is the 
only one : suppose that Xq is a minimum of J. Then J{xo) = and thus 

VAeM'' {\,xo)-L{X)<0 

Hence for alH > and A e M"^, 

m) - m ^ mi ^ {tx,xo} ^ 
t t - t ^ ' °' 

o 

Since L is differentiable at e Dl, letting t go to 0, we get 
VAeK'' (VL(0),A) > (a;o,A) 

It follows that xo = VL(0) = m. □ 

We notice that Cramer's theorem (see [5]) links J and the large deviations of 
{Xi + ■■■ + Xn)/n where {X„)neN is a sequence of real- valued independent 
and identically distributed random variables in R'^. This is why J is called the 
Cramer transform. 

A probability measure i/ on M is said to be degenerate if it is a Dirac point mass. 
We will generalize this definition for measures on W^. We refer to [4] and [11]. 

Definition 12. A probability measure v on , d>2, is said to be degenerate 
if its support is included in a hyperplane ofW^, i.e., there exists a hyperplane H 

ofW^ such that v{n) = 1. 

The following lemma illustrates the interest of this concept. 

Lemma 13. (a) Ifv is degenerate then its Cramer transform J vanishes outside 
of a hyperplane containing its support. 

(b) If Z is a random variable whose distribution is v, which is non-degenerate, 
then its covariance matrix Gz is invertible. 

Proof, (a) We assume that "H is the hyperplane given by (ao, z) = t, with t gM. 
and ao € M'^\{0}. We set = aoi/ll'»o||- We notice that z e H if and only if 
z — to belongs to the orthogonal of oq. Thus for any x ^H, 

J{x)= sup ((A,a:-io) -In / e'-^'"'-*"^ du{z)) 

> sup ( (A, a; — to) " In / e^ diy{z) 
= sup {k{ao,x — to)) = +00 
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(b) Wc have Tz = E{Y*Y) with Y = Z - E{Z). The matrix Tz is symmetric 
and thus it is diagonalizable. To conclude that Tz is invertible, it remains to 
prove that is not an eigenvalue . Suppose that it is the case : there exists a 
vector X y^O such that Fzx =* (0, . . . , 0). Then 

E(||*ra;||2) = ECxY*Vx) =*xE{Y*V)x =*xTzx = 

Therefore ||*i^a;||^ = almost surely and thus 

d 

XiYi = a.s 

1=1 

That is, with probability 1, 

Z & {z eM.'^ : {x,z) =E{{x,Z))} 

This is absurd since u is non-degenerate. Hence Tz is invertible. □ 

From now onwards, wc assume that is a non-degenerate probabihty measure 
in R''-. We are interested in the points A realizing the supremum defining J{x), 
for X e Dj. We denote by C the closed convex hull of the support of v. 

Lemma 14. Let u he a non- degenerate probability measure in R'^. The interior 
of C is not empty and C C Dj C C. Moreover for any x € C, the supremum 
defining J{x) is realized for some value \{x) G D^. 

Proof. The non-degeneracy of v means that its support is not included in a 
hyperplane of W^. Therefore the support of v contains d linearly independent 
vectors and the interior of the convex hull of these vectors is non-empty. Thus 
C is non-empty. 

We prove next the second assertion. We first show that Dj C C (see corollary 
12.8 of [5]). Suppose that C (otherwise the result is immediate). Let x ^C. 
By the Hahn-Banach theorem, there exists A € M'' and a e IR such that 

VyeC (A,y) <a< (A,a;) 

Since u{C) = 1, for any f > 0, 

J{x) > {tX,x) — In / exp{{tX, y}) di>{y) 

= — In / exp(t(A, y) — t{X, x)) dv{y) > t{{X, x) — a) 
Jc 

o 

Sending t to -hoc, we conclude that J{x) = +oo. Thus Dj c C. Let x & C and 
let (A„)„gN be a sequence in R'' such that 

J{x) = lim I (A„,a;) - In / exp((A„, z)) dv{z) ) 



= — In lim / exp{{Xn,z — x))di'{z) 

We suppose that |A„| — >■ +00 and we show that it leads to a contradiction. For 
all n € N, we set u„ = A„|A„|~"'^. Then (un)neN is a bounded sequence. Thus, up 
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to the extraction of a subsequence, wc might assume that it converges to some 
vector u G M'' whose norm is 1. Let v belong to the support of p and let U be 
an open subset of M"^ containing v. We have then i'(C/) > 0. Suppose that for 
any z G U, {u,z — x) > 0. Then, by Fatou's lemma, 

+00= / liminf exp( I A„ I z — a;)) 
Jb "-^+°° 

< liminf / exjp{\Xn\{un, z — x)) di'{z) 

n-)-+oo 

Hence 

exp(— J(x)) = lim / exp{\Xn\{un, z — x)) di'{z) = +00 

Thus J{x) = —00, which is absurd since J is a non-negative function. We 
conclude that for all v in the support of u and for any open subset U of M** 
containing v, there exists z G U such that {u, z — x) < 0. It follows that, for 
any v in the support of 1/, {u,v) < This inequality is stable by convex 

combinations, thus 

Vy e C {u, y) < {u, x) 




Case where v is discrete and charges five points of M.^ 

o 

Since x G C, there exists a ball B^; centered at x and contained in C. Thus there 
exists yo G such that {u, yo) > {u, x), which is absurd. Therefore (A„)„£n is a 
bounded sequence. Hence there exists a subsequence (A(^(„))„£n and X{x) G M** 
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< (A(x),x) 



such that A^(„) X{x). By Fatou's lemma, 

J{x) = {X{x),x) — hi lim / exp{{Xn, z)) dv{z) 

< (A(a;), a;) — In liminf / exp({Xn,z))dv'(z) 
In / liminf exp((A„, z)) dz^(z) 

= (A(a;), a;) — In / exp{{X{x),z))di'{z) 

< J{x) 

Thus J{x) = {X{x),x) — L{X{x)). Since L{X{x)) ^ — oo, this formula implies 

that J{x) < +00 and thus that C C Dj. Moreover if L{X{x)) = +00 then 
J{x) = —00, which is absurd. Therefore L{X{x)) < 00. This shows that the 

supremum defining J{x) is realized at a point X{x) such that A(A(a;)) < +00. □ 

, 00 
If Dl is an open subset of K then for^all {x,y) G Dj = C, the supremum 

defining J{x) is realized at some A(a:) € Dl- This is the case when the support 

of is bounded, and also for the distribution Vp when p is the Gaussian A/'(0, cr^), 

where we have then Dl = M.x] - 00, l/(2cr^)[. 

Now we study the smoothness of J. 

Notation. If f is a differentiable function on an open subset U o/R'^, we denote 
by Dxf the differential of f at x £U. If f is real-valued, we denote : 
★ D^/ its second differential at x G U (considered as a matrix of size dx d) 
■k V/ the function U — > R'^ such that 

\/x&U VyeK<^ (V/(ar),y)=D,/(y) 

We define the admissible domain of J : 

Definition 15. Let v be a non- degenerate probability measure on such that 
the interior of Dl is non-empty. The admissible domain of J is the set Aj = 
yL{DL). 

The following proposition states that Aj, the admissible domain of J, is an open 
subset of R'^, and that J is C°° on Aj. 

Proposition 16. Let u be a non-degenerate probability measure on M*^ such 
that the interior of Dl is non-empty. Let Aj be the admissible domain of J. 

o 

(a) The function VL is a -diffeomorphism, from Dl to Aj. Moreover 

Aj C Dj = {x eR"^ : J{x) < +00 } 

(h) Denote by X the inverse C°° -diffeomorphism ofVL. Then the map J is C°° 
on A J and for any x G Aj, 

J{x) = {x,X{x)) - L{X{x)) 

V J(x) = {VL)-\x) = X{x) and J = (d^ (^^l) 

J 00 
(c) If Dl is an open subset of Mr then Aj = Dj = C where C denotes the 

convex hull of the support of v. 
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The points (a) and (6) are proved in [1] and [4]. For the sake of completeness, 

wc reproduce the proof below. 

o o 

Proof, (a) We iinow that the function L is C°° on and that for any X £ Dl 
and i,j e {!,..., rf}, 

dXi^ expL(A) igd 

d'L ^ /m. ZiZje(>^'^^ dv{z) _ (4, ^,e<^--) dv{z)) z,-e<^--) dv{z)) 
dXiXj expL(A) (expL(A))' 

= ¥.{Zx,iZx,j) -nZx,^)nZxj) 

where Zx = \Zx,i, ■ ■ ■ , Zx,d) is a random vector in R'^ with distribution /j.^, 
whose density is 

z I — > exp( (A, z) — L{X) ) 

with respect to ly. Thus T>\L is the covariancc matrix of fi^. Moreover /i,^ has the 
same support as v and thus it is non-degenerate. Therefore lemma 13 implies 
that D^i is invertible. Hence, for all A € Dl, D^i is a symmetric positive 
definite matrix. It follows that, for any x gR'^, the equation 

VL(A) = X 

o 

has at most one solution A(.r) G Dl. Indeed, if there exist two different vectors 
Ai and A2 in Dl such that V-L(Ai) = V-L(A2), then the function 

V' : i ^ (VL(iAi + (1 - i)A2), Ai - A2) 

is and real-valued on [0, 1] and verifies V'(O) — ^(1) = 0- Rolle's theorem 
implies that there exists Iq € ]0, 1[ such that tp'{to) = 0, i.e., 

(VL(Ui + (1 - t)A2)) I , Ai - A2) = 

\at n=ta I 

Setting Ao = t^Xx -|- (1 - io)A2 G Dl and = Ai - A2 G K'^\{0}, we have 

{Yyx,L{v\v) = ^ 

This contradicts the fact that D^L is positive definite. Hence VL is a bijection 
from Dl to Aj with inverse function 

A : a; G 1 — \ X{x) 

O ° / \ 9 

Moreover L is C°° on Dl and for any X G Dl, Da(VL) = D^i is an isomorphism. 
Thus the inverse function theorem implies that V L is a C°°-diffeomorphism from 
Dl to Aj. 

(b) For x G Aj, we define 

: A G R"^ ^ {x, A) - L(A) 

o 

The map fx is differentiable on Dl and 

VA G Dl V/x(A) = X - VL(A) 
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Wc have shown in (a) that, for aU x E Aj, V/j;(A) = if and only if A = A(a;). 
Since fx is concave, its suprcnium is realized at A(x), that is 

J{x) = /.(A(x)) = (x, A(x)) - L(A(x)) 

It follows that A J C Dj and that J is C°° on Aj. Finally for any x G Aj, 
u e W', 

{VJix),u) = (u,A(x)) + {DxX{u),x) - (VL(A(x)),D^A(u)) = {u,X{x)) 
since VL(A(a;)) = x. Hence 

Vx e Aj VJ{x) = X{x) = (VL)-^(2;) 
Differentiating VL(A(x)) — x, we get 

DA(.)i°D'^ = DA(x)i°D,A = Id 
whence the expression of J since D^^^^L is an isomorphism. 

, o o 

(c) If Dl is an open subset of R then lemma 14 implies that for a; e C = Dj, 
the supremum defining J{x) is realized at some point X{x) £ -Dl = Dj,. The 
function L is differentiable at A (a;) and the point (6) yields that 

X = VL{X{x)) e A{Dl) = Aj 

o o o 

Thus Dj c Aj. Finally, since Aj c Dj and Aj is open, we have Aj = Dj = C. 
This proves (c). □ 

Let V be a probability distribution on having a density with respect to the 
Lebesgue measure and let {Xn)n>i be a sequence of independent and identically 
distributed random variables with distribution v. The following theorem states 
that, under some hypothesis allowing the Fourier inversion, the density of the 
distribution of {Xi H + Xn)/n is asymptotically a function of 



J : a; e ^ 



sup ({t,x)-\n I e<*'^> ) 



We propose a proof, extracted from the article of C. Andriani and P. Baldi [1]. 
It relies on proposition 16. 

Theorem 17. Let v he a non-degenerate probability measure on M''.^ We denote 
by L its Log-Laplace and by J its Cramer transform. Suppose that Dl ^ and 
that there exists no > 1 such that 



V 



G Li(M'') 



We denote by Aj the admissible domain of J. For any x € Aj, we set jix the 
probability measure on such that 

exp{y + X, X{x)) , . , ^ 
expL(A(a;)) 

(where X is the inverse function ofVL). For n large enough, the Fourier trans- 
form of jjLx belongs to L"(M''). Let {Xn)n>i be a sequence of independent and 
identically distributed random variables with distribution v. For any n>no, the 
random variable Xn = (^i + • • • + Xn)/n has a density Qn with respect to the 
Lebesgue measure on satisfying : 
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(a) For X G Aj and for n large enough, 

(b) If Kj is a compact subset of Aj then, uniformly over x G Kj, when n goes 
to +00, 



The proof requires some preliminary results, which are presented next. 

Lemma 18 (Uniform dominated convergence theorem). Let X be a separable 
space and let {fl,J^,^) be a measurable space. Let f and /„, n > 1, be real or 

complex-valued measurable functions defined on X x Q. Suppose thai, for any 
u) € fl, the functions x i — > f{x,oj) and x i — > fn{x,u)), n G N, are continuous 
on X and that 

sup |/„(a;,a;) - /(a;,a;)| — > 

Suppose also that there exists a non-negative and integrable function g on Cl 
such that 

VneN yxeX VweO |/„(x,a;)| < g(a;) 
Then for any x G X , the function ui i — > f{x,ui) is integrable and 



sup 



/ fn{x,uj)diJ,{uj) - / f{x,ui)diJ,{co) 
Jn Jn 



— > 

n— >oo 



Proof. We adapt the proof of the classical dominated convergence theorem 
in [16]. Sending n to +oo in the domination inequality, we get 

\/{x,u!) € X X Q. \f{x,u!)\<g{u!) 

This shows that oj i — > f{x, co) is integrable. For any n € N, we set 

hn-.OJl > sup \fn{x,U)) - f{x,U))\ 

For all n G N and u E fl, the function x E X i — > \fn{x,uj) — f{x,u!)\ is 
continuous and, since X is separable, its supremum is equal to its supremum 
on a countable dense subset of X. Therefore /i„ is a measurable function. We 
have that {2g — /i„)„gN is a sequence of non-negative functions whose limit is 
the function 2g. Fatou's lemma implies that 

/ 2gdii = / liminf {2g — hn) dfi < liminf / {2g — hn) d^i 

Jn Jn "^+°° n^+oo 

2g dji — limsup / /i„ dji 

rn-l-oo Jn 

Since g is integrable, we get that 

limsup / hnd/j, <0 

n^+oo Jn 
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Hence hndfj, ^ since for any n e N, /i„ is a non-negative function. Finally 



sup 

xex 



fn{x,uj)dti{uj) - / f{x,uj)dii{uj) 



< sup / \fn{x,Uj) - f{x,Uj)\dll{oj) 



< 

/a 



/ hndji - 
Jo. 



and the lemma is proved. □ 

Lemma 19. Let v he a probability measure on M''. We denote by L its Log- 
Laplace. Let K be a compact subset of Dl ■ Then the function 

{s,t)^ M{s + it)= [ e^'+'^-^'Ui^ix) 
is uniformly continuous on K xR'^. 

Proof. For n > 1, we denote by B„ the open ball of radius n centered at the 
origin and we set 

fu-.y^ I e<«'^> dv{x) 



The sequence (/„)«>! is a non-increasing sequence of continuous functions on 
K C Dl , which converges to the null function. By Dini's theorem, the sequence 
(/ra)n>i converges uniformly to the null function on K. Let e > 0. There exists 
no > 1 such that 



We define next 

yx,s,teR^ g{x,s,t) = exp{{s + it,x)) 

The function g is uniformly continuous on B„(, x K x (its differential is 
bounded on this set, hence g is lipschitz). Thus there exists S > such that for 

s,ue K and t,v e K"*, 



{s,t) - {u,v)\\ < S =4- VxeB„o \g{x,s,t) - g{x,u,v)\ < - 



£ 

2 

Therefore 



\M{s + it) - M{u + iv)\ 

< / \g{x,s,t) - g{x,u,v)\di>{x) +2 sup / e^''^''^^ dv{x) 

This proves the uniform continuity of (s, t) i — > M{s + it) on K x M.'^. □ 

We will use the Riesz-Thorin theorem to prove our last lemma. Recall that the 
norm of a continuous linear operator T from L^'(M'^) to L'(R'^), with {p,q) € 
[1, +oo]^, is defined by 
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Theorem 20 (Riesz-Thorin). Letpo, pi, qo andqi in [l,+oo] such thatpo ^ pi 
and Qo ^ qi- For any t e [0, 1], we put 

1 l-t t , 1 1-t t 

— = 1 ana — = 1 

Pt Pq Pi qt qo qi 

Let T be a continuous linear operator from LP°(R'^) to L'^°{R'^) with operator 
norm Mq. Suppose that T is also continuous from Lp^{U.'^) to L'^^{U.'^) with 
operator norm Mi. Then, for any t G [0,1], T is a continuous linear operator 
from LP* {R'^) toLi*{R'^) and 

V/ e L^'*(M<^) ||T(/)||,, < Mt*Ml\\f\\,, 

We refer to chapter 1 of [3] for the proof of this theorem. We apply next the 
Riesz-Thorin theorem to the Fourier transform : 

The map which associates each intcgrablc function to its Foiirier transform is 
a continuous hnear operator from L^(M'') to L^(M'') with operator norm 1. 
Moreover Plancherel theorem guarantees that 

V/eL2(IR'^) ll/lb = (27r)'^/2||/||2 

For any t G [0, 1], we put pt = 2/(2 — and qt = 2/t. We easily check that pt 
and qt are conjugate and that 

1 l-t t , 1 l-t t 

— = 1 and — = 1 

Pt Po Pi qt qo qi 

with Po = 1, pi = 2, qo = +oo and qi = 2. The Riesz-Thorin theorem implies 
that 

V/ e LP*(E^) 11/11,, < li-*(27r)"/2||/||p, = (27r)'^/«'||/||p, 
Since pt E [1,2], we get the following inequality : 

Lemma 21 (Hausdorff- Young inequality). Let d> 1 and p £]1,2]. We denote 
q = p/{p - 1) e [2, -|-oo[. Then 

V/GL^'(M'^) 11/ II, < (27r)'^/'||/||, 

Proof of theorem 17. We denote by ip the Fourier transform of v. We have 

Vn > 1 = z?" = y," 

By hypothesis, </9"° e Li(M'^) and 

Vn > no |i^| = Iv-I"" |(^|"-"° < |(^r° 

Thus V*" e Li(M'^) and the Fourier inversion formula implies that i/*" has a 
density /„ given by 

Vx G M fn{x) = [ e-^<*'-V" W dt 
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We also get that /„ e L°°(]R''). Let s £ Dl and n > 1. The function x i — > 
^^^'^^ fn{x) is non-negative and its integral over R'^ is M(s)" < +00. J^et us 
denote by (fis,n its Fourier transform. Let p €]1,2] be such that ps G Dl- We 
have 

which is finite. The Hausdorff- Young inequality implies that (ps,no G L'^(R'') 
where q = p/{p — I) € [2, -|-oo[. However 

e R'^ ^s,n{t) = I e^^*'^>e<"'^V„(a;) dx = M{s + it^ 

Thus for any n > qno, 

[ \^s,n{t)\dt= I |M(s + it)r<'9|M(s + ii)|"-"'"'di 

< M(s)"-"°« / (|M(s + dt 

= M(S)"-""'' / \^,^no\^ dt < +00 

Thus (ps,n G L^(M'') and the Fourier inversion theorem yields that for all a; G R, 

= 7^ [ dt=-^J e-^<*'->M(5 + ii)" 

Moreover for x S and t e R*^, 

c!ny + a;j-e M(A(a;)) 

where we made the change of variables 2; = ?/ + a;. It follows that € L" 
for n > nog. Notice that 

Va; e R"^ ff„(a;) = n'^fn{nx) 

Therefore 

If a; e j4j then proposition 16 implies 

J(a;) = (A(a;),a;)-L(A(a;)) 
thus, applying the above inequality to s = A(a;), we get 

^(^Y f (e-<A(x)+i*,x>+L(A(x)+it))" 

V27r/ J^d ^ ^ 
= ^iL^'^e""'''''^ y" (^e-^i=>')-iKx)+it,x)+L{X{x)+it)y 

^ (^Jl.y f,~nJ{x) f (^^-i{t,x)-L{\(x))+L{X(x)+it)y 
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This equality is valid when n > noq. This proves the point (a). 

Now let us prove the point (b). Let Kj be a compact subset of Aj. We notice that 
q depends on x G Aj, but, by compactness of Kj, we can choose q uniformly 
over X e Kj. For any x G Kj, the mean of /U^ is 

/ y , , diy(y + x) = / (z - x)—-—-- di^(z) =VL(X(x)) - X = 

and its covariance matrix is = D^^^^L since for 1 < i,j < d and s G Dl, 

/r. yi%e<^(-)'^+-> du{y + x) _ J^,{zi - Xi){zj - x,)e<^(-)'-> du{z) 



M{X{x)) M{X{x)) 



M{X{x)) ' ^ dsiSj 

When t ^ 0, uniformly over x G Kj, 

1 



(A(x)) 



l^{t) = 1 - -{T^t,t) + o{\\t\\^) 
Indeed 

(x,t) ^ fi,[t) ^ e M(A(x)) 

is C°° on A J X R'' (by proposition 16), thus the Taylor-Lagrange formula gua- 
rantees that the remainder term is uniformly controlled over x G Kj. Therefore, 
for any t e IR"^, uniformly over x G Kj, 



l^x\A=\ — ^ exp 



n—^oo 



-^(rxt,t)) 



The functions x i — Jl^ and x i — cxp {—(Txt^t) /2), t E M."^, are continuous on 
Kj. In order to apply the dominated convergence theorem (the uniform variant), 
we need to get a uniform domination of the sequence of functions. For x G Aj, 
Tx is a positive definite symmetric matrix thus Sx, its smallest eigenvalue, is 
positive. The largest eigenvalue of the inverse of is e~^. Therefore, for any 
X G Aj, 



Ex = (max { a : a eigenvalue of ^ }) = ( sup — ^ ' ^ | 



-1/2 



The term on the right is the inverse of the operator norm of the linear application 
associated to the matrix F"'^. Moreover x i — > Yx — D^^^ji is continuous on Aj 
thus the function x i — > Sx is continuous. Let us denote by £o its minimum 
on Kj. The compactness of Kj enstires that eo > 0. The previous expansion 
implies that there exists ^ > such that 

V(t, x) G B(0, 6) X Kj |fc(t)| < 1 - 1 ^ I^F, - ^Id) t, t 

The spectral theorem for real symmetric matrices yields that, for any x G Kj, 
the matrix F^, — eolrf is positive symmetric. Thus 

Vt G K-^ ((f, - |ld) t,t) - = ((F, - Sold)t,t) > 
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It follows that 



y{t,x) e B{0, S) X Kj \lli{t)\ < 1 - ^\\t\\ 



Since I — y < e ^ for all y > 0, we get 
Vn > 1 V(i, x) e B(0, 6y/n) x Kj 



The right term is integrable and does not depend on a; e Kj and n. The uniform 
dominated convergence theorem (lemma 18) implies that, uniformly in Kj, 



i 



(J-V 



t\\<S^ ^V"' 



dt ^ I exp(-i(r^i,0) dt 



Moreover this second integral is equal to (27r)''/^ (detr^;) ' and proposition 16 
guarantees that for x E Aj, D^^^-^L is the inverse matrix of J. Hence, when 



n — 7- oo, uniformly over x G Kj, 



[ iTSTdt = 71-^^/2 f fty ^ f^Y' (detD^j) 

J\\t\\<s A\t\\<sV^ ^^/nJ \nJ ^ 



1/2 



Let us focus on the remainder of the integral. Wc set 

h: X G Kj I — > sup 

ll«ll>'5 

The function A is continuous thus X{Kj) is compact and lemma 19 states that 
the function (s, t) i — > M{s + it) is uniformly continuous on X{Kj) x W^. The- 
refore the function 

XI — > sup \M{X{x) +it)\ 

m>s 

is continuous on Kj. However 



Va; e Kj h{x) = sup |/Xx(i)l = 

u\\>s 



— -— - sup \M{X{x)+it)\ 
M{X{x)) ut\\>S 



Hence h is continuous on Kj. By compactness of Kj, there exists xq G Kj such 
that 

sup sup \]j^{t)\ = sup h{x) = h{xo) = sup \]j^{t)\ 
xeKj\\t\\>s xeKj \\t\\>s 

Finally, just like u, the law has a density and the Riemann-Lebesgue lemma 
implies that 



\\t\\^+o 







Moreover lemma 4 of chapter XV. 1 of [11] guarantees that for any t ^ 0, 
< 1. Therefore there exists k € ]0, 1[ such that 



sup 

l|t||>5 



We get 



sup sup |Atx(i)l < K < 1 

xeKj \\t\\>s 
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It follows that for any x € Kj and n > noq, uniformly over x G Kj, 

I mrdt<l \fi;{t)rdt<K"-'''"^ [ \ii;,{t)r'dt = o{i) 

J\\t\\>5 J\\t\\>5 J\\t\\>6 

since k <1 and Ji^e L"«9(R''). Finally 

9n{x) = ( [ rAtr dt+ [ KAtr dt] 

^^TT/ \"'||t||><5 "'Pll<<5 / 



= (^)'"%-""'(<i«Dy)"'(l+»(l)) 



The expansion is uniform over x G Kj^ by the previous results and the boun- 

1 /2 

dedness of x i — > (det j) on Kj, since it is continuous. This ends the proof 
of theorem 17. □ 

Proposition 22. Let v he a non-degenerate probability measure on Mf^ such 
that Dl 7^ 0- If there exists m e N and p €]1,2] such that v*"^ has a density 
fm € LP(IR'^) then the hypothesis of theorem 17 are verified. 



Proof. It follows from the Hausdorff- Young inequality that fm S L''(IR"), with 
f = p/{p ^ !)• Moreover is bounded thus fm S L5(IR'^), where g is a positive 
integer larger than r. Therefore 

Hence the hypothesis of the theorem are verified with no = mq. □ 

7 Minima of I-F 

Let p be a probability measure on K. We define 

A : {u,v) e ^ In [ 6"^+"^' dp{z) 
Js. 

We denote by I the Fenchel-Legendre transform of A. 

In this section, we consider the minima of the function I—F when p is symmetric. 

a) Admissible domain of / 

We begin by giving some properties of I, which arc consequences of the results 
stated in the previous section. Let be the distribution of {Z, Z"^) when Z is 
a random variable with law p. We suppose that the support of p contains at 
least three points so that z/p is a non-dcgcncratc measure on M^. The function A 
is the Log-Laplace of and its domain of definition contains Mx] — oo, 0[, 
thus its interior is non-empty. The function I is the Cramer transform of Vp and 
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o 

we denote by Aj = VA(Z)a) the admissible domain of /. Proposition 16 implies 

that : 

o 

(a) The function VA is a C°"-difFeomorphism from Dj^ to Aj. Moreover 

Ai c Di = {x : I{x) < +00 } 

(b) The fimction / is C°° on Ai. If i — > {u{x,y),v{x,y)) is the inverse 
function of VA then, for any (.t, y) G Aj, 

I{x, y) = xu{x, y) + yv{x, y) - A{u{x, y), v{x, y)) 

VI{x,y) = (VA)-^(a;,y) = {u{x,y),v{x,y)) 

9 O O 

(c) If Da is an open subset of R then Aj = Dj =C where C is the convex hull 
of the set { (a;, a;^) : x is in the support of p }. 

b) Minimum of / — F on A* 

Let p be a symmetric and non-degenerate probability measure on M. Jensen's 
inequality gives us 

V(«, v) e In / 6""+""' dp{z) > I {uz + vz^) dp{z) = va^ 

thus 7(0, cr^) < 0. Since / is non-negative, then 7(0, cr^) = and 

inf(7-F)e[-^,0] 
The function I is even in the first variable. Indeed, if (a;, y) G M^, then 
I{-x, y) = sup ( -xu + - In / e"''+'"'^ dp{z) ) 

= sup ixu + yv -hi I 6"""'+'"'^ dp{z) ) = I{x, y) 

Assume that 1 — F has a unique minimum (a;o, j/o) on ^ ■ Then (— a^o, 2/o) is also 
a minimum oil — F since 

-''(-a;o, 2/o) - F{-xq, yo) = I{xo,yo) - F{xo, yo) 
The uniqueness of the minimum implies that a;o = so that 

inf {I-F) = 7(0, yo) - F{0, yo) = 7(0, yo) > 

Since 7(0, ct^) = we have yo = a^. 

In this section, we will show that, if p is symmetric, then I — F has a unique 
minimum on A*, which is at (0,(7^). 
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Consider first the case of a Bernoulli distribution. Let c > 0. Suppose that 
p = ((5_c + <5c)/2. The law p is centered and its variance is . For all {u,v) G M^, 

K{u, = In f 6"^+"^' dp{x) = vc^ + Incosh(Mc) 

For any {x,y) e M^, by studying the function {u,v) i — > xu + yv — A{u,v), we 
can determinate its supremum. We get that / is finite on Dj = [— c, c] x {c^}, 
/(— c, c^) = /(c, c^) = In 2 and for any x g] — c,c[, 

I{x, c^) = -z- {{c + x) ln(c + x) + {c — x) ln(c — a;)) — Inc 

The function g : x i — > I{x, c^) - x^/{2c'^) is on ] - c, c[ and 

Vx e ] — c, c[ ^'(ic) = - ^arctanh ^— ^ — — ^ 

c — X'^ c 

Thus g' is non-decreasing on [— c, c] and, since ff'(O) = 0, it follows that g has 
a unique minimum at 0. Therefore I — F has a unique minimum in A* at 
(0,c2) = (0,a2). 

The previous results yield the following lemma : 
Lemma 23. Let c > 0. We define 



The function 



ipc ■ X gR I — > sup ( ua; — In cosh(uc) 

mSR 



x G [-C, c] I — > ipc{x) - ^ 



has a unique minimum at and ific{0) = 0. 

Notice that the Bernoulli case is special since, if X is a random variable with 
distribution p = (5-c + <5c)/2, then X"^ = (? almost surely. Thus 



-Z'^ U x\ + ---+xl ) ^{^l+-+<>o} 11 "^P^^^^ 



1 / {Xi +--- + Xn)^\ TT J / ^ 

i=l 

This is exactly the classical Curie- Weiss model. 

In the following, we suppose that the support of // contains at least three distinct 
points. We first show that, if Da is an open subset of M^, then I—F has a unique 
minimum at (0,0-^). 

In the subsection a), we saw that, if the support of v contains at least three 
distinct points and Da is an open subset of M?, then / is diflerentiable on the 
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interior of its domain of definition Di and, if {x,y) i — > {u{x,y),v{x,y)) is the 
inverse function of VA, then 

o dl 
y{x,y)GDi —{x,y)=u{x,y) 

If we show that u{x,y) > x/y for a; > and y > 0, then, by integrating this 
inequahty and using the fact that I is even in its first variable, we get that 

2 

V(a;, y) G Dj I{x, y) - J(0, y)>^= F{x, y) 

To obtain that I — F has a unique minimum at (0,cr^), it is enough to extend 
this inequality to the boundary points of Dj (if they exist) and to show that 
/(O, .) has a unique minimum at cr^. 

The following lemma is the key result to establish the uniqueness of the minimum 
oi I — F, when p is symmetric. 

Lemma 24. Let p be a symmetric probability measure whose support contains 
at least three points. We have u{x, y) = if x = and 

X 

u{x,y) > — if X > 

X 

u{x, y) < — if X <0 

Proof. The vector {u,v) = {u{x,y),v{x,y)) verifies 

[X, y) - VA(«, v)-^j^ ^„,+,,. , ^^^^^ 

The distribution p is symmetric, thus 

/• r+oo 

/ 26""+""' dp{z) = / 2zsmh{uz)e'''' dp{z) 
Jr Jo 

This formula shows that u and x have the same sign. Moreover for any z > 0, 
tanh(z) < z thus, if a; > then sinh(u^;) < u2;cosh('U^;). Therefore, using the 
symmetry of p, 

/+°° 2z\oshiuz)e"'' dp{z) L z2e«^+'-' dp{z) 
X < u — T ■ — , , , ^ = u ■ — ,j , , , = uy 

J^^uz+vz^ dpiz) J^e--+--' dp{z) 

Since a; > 0, u > and y > 0, we conclude that u > x/y. Similarly, we show 
that if a; < then u < x/y. □ 

We can now prove the following inequality : 

Proposition 25. If p is a symmetric probability measure on M with variance 
cr^ > and such that D\ is an open subset of M.^ then 

\f{x, y)GRx M\{0} I{x, y)-^> 1(0, y) 
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Notice that this rcsuh encompasses the case of a symmetric measure with boun- 
ded support, because in this case Da = K^. In proposition 31, we shall extend 
the inequality to any symmetric distribution on M. 

Proof. We already treated the Bernoulli case (see lemma 23). We assume next 

that the support of p contains at least three points. Wc denote by Ai the ad- 
missible domain of / and by C the convex hull of the set 

{ (x, x^) : a; is in the support of p } 

In the subsection a), we saw that, if Da is an open subset of R^, then Ai = 
Dj = C C A*. Moreover / is C°° on Dj and if {x,y) i — > {u{x,y),v{x,y)) is the 
inverse function of VA then 

o dl 
V(x, y) eDi —{x,y)=u{x,y) 

Let us examine the structure of the set Dj. We put 

\/y>0 Dj^y = n (M X {y}) 





A 






Di 


= cJ 


\ ' 


y 

(.-- 


-Jhi 



y = x^ 





Case where p is symmetric discrete and charges 5 points 

o 

Let y > be such that {x, y) G Dj for some a; S M. The set Di^y is a convex 
subset of M. Moreover x \ — > I{x, y) is even, therefore Di^y is an open interval 
] - a{y), a{y)[ with a{y) G [0, We have 



\/xeDi,y Iix,y) - l{0,y) 



u{t, y) dt 



Lemma 24 implies that for any t > 0, u{t,y) > t/y. By integrating and using 
the fact that / is even, we get that 

2 

VxGDj^y I{x,y)-I{0,y)>^ 

o 

and there is no problem of definition at y = since Dj c A* does not contain 
M X {0} and Difl = 0. Moreover 



I{x,y)-IiO,y) 
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is non-dccrcasing on Z?/.j,\{0} since / is convex. Therefore, if ~a{y) and a{y) 
belong to Dj^y, then the previous inequahty extends to a; = —a{y) and x = a{y). 
We have shown that 



X2 



•^{x,y)&Di V2/>0 I{x,y)~I{Q,y)> — 

except for the points {x,y) of the superior and inferior borders of if they 
exist. More precisely, we set 

= inf { a;^ : x is in the support of p} > 

and 

= sup {x^ : X is in the support of p} < +00 

If if = and L = +00 then the inequality is already proved on £)/\{(0, 0)}. 
Suppose that K'^ > 0. Let y = K'^ and a; e M. We define 

/ : (u, w) e I — >ux + vK'^ - A{u, v) 
Denoting ck = p{{K}), we have for all {u,v) G M^, 

f{u, v)=ux- ln{2cKC0sh{uK)) -In [ ^uz+v{z^-k^) ^^^^^j 

For any z G K\] — K, K[, the function v 1 — > exp{v{z^ — K^)) is non-decreasing. 
Therefore 

sup f{u, v)=ux- ln{2cKC0sh{uK)) - In ( lim / ^uz+v{z^-k^) ^^^^-j j 

iJGR \^^~°°^R\[-is:,is:] / 

= ux — ln(2cKC0sh(uii')) 

by the dominated convergence theorem. Indeed 



yzem.\[-K,K] yv<-i 



<- ^uz-iz^-K^) 



and the map z G M\[—K,K] 1 — > e"^ ^ Ms integrable with respect to p 
since it is bounded (it is continuous and goes to when \z\ goes to +00). Hence 

I{x, K^) = sup f{u, v) = sup {ux — ln(2ci<-cosh(u/"ir)) } 

In fact, we come back to the Bernoulli case. The reason is that, if we condition 
on Tn = in our model, then for any i, X^^ — ~K or K . 

li Ck > 0, then lemma 23 implies that for all x & {—K,K\, 



I{x,K^)-I{0,K^) = ipK{x)> 



,.2 



2ii'2 



If Ck = then for any x ^ 0, I{x, K^) = +00 so that the inequality is verified 
for y = K'^. 

If L < +00 then we show similarly that for all x G [—L, L], 



7(x,L2)-/(0,L2)>^ 
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Therefore for any {x,y) G £)/\{(0,0)}, 

I{0,y)<I{x,y)-^ 

Notice that for any j/ € R, by the convexity and the symmetry of x i — > I{x, y), 
if l{0,y) = +00 then for all x ^ 0, I{x,y) = -\-qg. Therefore the inequality 
extends to each subset of which does not contain M x {0}. □ 

In the previous proof, if we take x = y = 0, then for any m e R, the function 
V I — > K{u, v) is non-decreasing on R. Therefore 

infA(u,v)= lim A{u,v) = lim ( lnp({0}) + In / e"^+"^' dp(^) ) 
vm y Jr\{o} ) 

= lnp({0}) 

by the dominated convergence theorem. Hence 

inf AKz;) = inf (lnp({0})) = lnp({0}) 

This is valid for any probability measure p in R. This yields the following lemma : 

Lemma 26. If p is a probability measure on R then /(0,0) = — lnp({0}). 

If Da is an open subset of R^ then (0, 0) G Da = Da- It follows from the point 
(d) of proposition 11 that /(O, .) has a unique minimum at cr^. Therefore, the 
inequality of proposition 25 implies the following corollary : 

Corollary 27. If p is a symmetric probability measure on R with variance 
cr^ > and such that Da is an open subset of M? then the function 

ix,y)eA*^I{x,y)-xy{2y) 

has a unique minimum at {0,a^). 

Now we will extend this result to any symmetric probability measure such that 
(0, 0) e Da- For this we need Mosco's theorem, which we restate next. 

Definition 28. Let f and fn, n gN, be convex functions from M.'^ to [—00, +00]. 
The sequence (/n)neN is said to Mosco converge to f if for any x G R"^, 
■k for each sequence {xn)neN ^'^ converging to x, 

liminf /„(a;„) > f{x) 

* there exists a sequence {xn)n&i in R'' converging to x and such that 

limsup/„(a;„) < f(x) 

n— >+oo 

We write then 

f f 

In ' / 

n— >oo 
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If / is a convex function from M.'^ to [—00, +00], we define its Fenchel-Legendre 
transform /* by 

yxem'^ r{x) = snp{{t,x)-f{t)) 

Theorem 29 (Mosco). Let f and fn, n G N, be convex functions from M'' to 
[—00, +00] which are convex and lower semi- continuous. We have the equivalence 

e p „* 

In ' / — ^ Jn ' / 
n— >oo n— >oo 

We refer to [13] for a proof. 

Proposition 30. Let v he a probability measure on . We denote by L its 
Log-Laplace. Let {K„)n£fq be a non- decreasing sequence of compact sets whose 
union is W^. For all n E N, we set Vn — v{.\Kn) the probability v conditioned by 
Kn and we denote by L„ its Log-Laplace. Then 

Ln — > L 

n— >oo 

Proof. For n large enough, the compact set Kn meets the support of and we 
have for any borel set A, 

Thus, for n large enough and A S K'', we have 

L„(A) = In / e<^'^> dvKn{z) = In / e<^'^> dv{z) - \nv{Kn) 

By the monotone convergence theorem, 

lim L„(A)=ln/ lim (1k„ (^)e<^'^>) - lim \nv{Kn) = L{X) 

Hence the second condition of Mosco convergence (with the limsup) is satisfied 
with the sequence {Xn)nm constant equal to A. 

Let A e K'' and (A„)„£n be any sequence converging to A. Fatou's lemma implies 
that 

expi(A)= / liminf lK„(-2)e<^'"^>rfz/(^) < liminf / lK„(^)e<^'"'=> ^1/(0) 

Jgd n-i.+oo J^i 



Therefore 



L{X) < liminf (-L„(A„) +lnz/(/r„)) = liminf L„(A„) 



Thus the first condition of Mosco convergence (with the liminf) is verified. □ 
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Proposition 31. Let p he a symmetric probability measure on M with variance 
cr^ > 0. We have 

V(a;,t/)eA* 7(ar,y)-^>/(0,y) 
Moreover, if A is finite in a neighbourhood of (0, 0), then the function 

{x, y) G A* I — > I{x, y)- 
has a unique minimum at (0,0-^) where it is equal to 0. 

Proof. For any n G N, we put Kn = [— n, n]^. For n large enough so that Kn 
meets the support of Vp, we define Un = T^p{-\Kn)i A„ its Log-Laplace and /„ its 
Fenchel-Legendre transform. For all {u,v) G M?, 



An{u, v)=\n [ e""+"* di^p{s, t) - \nVp{Kr. 



) <K{u,v)-\nyp{Kn) 



Applying the Fenchel-Legendre transformation, we get 

Vt/eM I{Q,y)<Ir,{G,y)-\nVp{Kr,) 

Moreover the measure z/„ has a bounded support thus proposition 25 and the 
previous inequality imply that 

-2 



Va;eM Vy>0 7(0,y) + ^ < 7„(a;,y) - ln!/p(Jr„) 



It follows from proposition 30 that (A„)„gN Mosco converges to A. Hence, by 

Mosco's theorem, {In)nen Mosco converges to I. In particular, for € 
such that 2/ > 0, there exists a sequence {xn,yn) G converging to {x,y) and 
such that 

limsup7„(a;„,y„) < I{x,y) 

n— >+cx) 

Since y > 0, there exists no > 1 such that t/n > for all n > no- Therefore 

x^ 

\/n > no l{0,yn) + TT- < In{Xn,yn) " lni'p(ii:„) 

^yn 

Moreover Vp{Kn) — ^ 1 when n ^ oo. Hence 

x^ 

limsup /(O, y„) + — < I{x, y) 

n—¥-\-oo 



Finally / is lower semi-continuous, thus 



liminf/(0,y„) > J(0,y) 

n— f+oo 



It follows that for {x, y) e such that y ^ 0, 

/(0,y) </(x,2/)- 



2y 



Suppose in addition that A is finite in a neighbourhood of (0,0). Point [d) of 
proposition 11 implies then that 7(0, .) has a unique minimum at cr^. Therefore 
{x,y) G A* I — > I{x,y) — x^ /{2y) has a unique minimum at (0, cr^) where its 
value is 0. □ 
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c) Expansion of I — F ciround its minimum 

If p is a symmetric probability measure whose support contains at least three 
points and if (0,0) e Dl then {0,a'^) = VA(0,0) e VA(L'a) = Ai, the admis- 
sible domain of /. We saw in subsection a) that / is C°° in the neighbourhood 
of (0,(7^) and that 

VI{0,a^) = {u{0,a%v{0,a^)) = (VA)-i(O, cr^) = (0,0) 

since D^q g^A is the covariance matrix of t'p (see the proof of proposition 16). 
Up to the second order, the expansion of / in the neighbourhood of (0,cr^) is 

x"^ {y - (T^)^ 
2^ ^ 2(/i4 - 0-4) 



The expansion of F up to the second order in the neighbourhood of (0, cr^) is 

Therefore, in the neighbourhood of (0,cr^), 

l{x, y) - F{x, y) = + o{\\x, y - a^f) 



We need to push further the expansion of I — F. 

Consider the case of the centered Gaussian distribution with variance cr^. We 
can compute explicitly / : 

y{x,y) e A* I{x,y) = 1(^1^-1 -In (^^)) 
If ip is the function defined in proposition 8, then 

I{x,y)-F{x,y)=i;(^,^)-l 
Lemma 9 implies then that, in the neighbourhood of (0,0"^), 

/(.,,)- F(.,,) ~ -(-, + - i) ) = ^ + 

In fact, we have a similar expansion in a more general case : 

Proposition 32. If p is a symmetric probability measure on M. whose support 
contains at least three points and such that (0,0) e D\ then I is C°° in the 
neighbourhood o/(0, cr^). If Ha denotes the fourth moment of p then, when {x,y) 
goes to (0, a^), 

x^ {y-a^f , p.Ax'^ 
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Whenever p = J\f{0, cr^), we have 

2(/i4 - (T*) = 2(3(7^ - a^) = 4(7* 

and 

M4/(12ct^) = 3aV(12c7^) = (4(7*)-^ 
This is what we obtained before the proposition in the Gaussian case. 
Proof. If (0, 0) e Da then 

(0, a^) = VA(0, 0) e VA{Da) = Aj 

The function / is on Ai and, if (a;,y) i — > {u{x,y),v{x,y)) is the inverse 
function of VA then, for all {x, y) G Aj, 

I{x, y) = xu{x, y) + yv{x, y) - k{u{x, y), v{x, y)) 

VI{x,y) = {VA)-\x,y) = {u{x,y),v{x,y)) 

o 

Moreover the hypothesis (0,0) G D\ implies that p has finite moments of all 
order. The expansion of F to the fourth order in the neighbourhood of (0, a^) is 

Therefore, in the neighbourhood of (0, 0), 

I{x, h + a'^)- F{x, h + a"^) = — j- + aa^x^ + 02,1x^/1 + ai,2a;/i^ + 

2(^4 - cr j 

+ a4fiX* + a^ ix^h + 02,2^;^^-^ + ai^^xh^ + 00^4/1^ + o(||a;, h\\^) 
with, for any e N such that i + j G {3, 4}, 

1 9*+^ J 2 

except for 

"^■^ = 2 ) + "^-^ = ia^w^"'" ^ " ^ 

If we prove that 04.0 > then the terms x/i^, It' , x^h, x^h^, xh^ and h* are 
negligible compared to 04,02;^ + ao,2^'^ when {x,h) goes to (0,0). Next, the 
symmetry of / — F in the first variable implies that 03,0 = 0. If we show that 
(i2,i = then we get 

I{x, y) - F{x, y) = (^^^—^^ + a4,oa:*) (1 + o(l)) 

when {x,y) — )• (0,C7^), so we have the desired expansion. 
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To conclude it is enough to show that 02,1 = and 04,0 = 1^4/ {12a^), that is 

(0, ^ and ^(0, = ^ 



For any j £ N, we introduce the function fj defined on Da by 



\/{u,v) € Da fj{u,v) 



Jr 



These functions are C°° on Da and they veriiy the following properties : 
★ /o is the identity function on M? 

For all j G N, fj{0, 0) = fij is the j-th moment of p. It is null if j is odd, since 
p is symmetric. Moreover 



, dA ^ , dA 

h = and 12 = 

ou ov 



★ For any j G 



^ = - fjh and ^ = /,+2 - /,/2 
Therefore, for all {x,y) G Aj, 

Denoting by 5 = {f2- fl){f 4- fi)- {h- f 1/2 f , the determinant of the positive 
definite symmetric matrix D^A, we get that for any {x,y) G Aj, 



)2 T ]_ ( 74-/2 /1/2 - fs 

g{u{x, y),v{x, y)) V/1/2 - /s /2 - /f 



Moreover (m(0, a^), t;(0, a^)) = (0,0) thus 

— 0,(7^ = 0,^2 = 0,0 = -i^ = — 

OX ox^ g a^{fM-(y ) <J 

-(0,<72) = — (0,a2) = ^^^-^(0,0) - 



dy ' dy'^ ' 5 ' a^ip^A - (t^) Pi - 

-(0,a ) = -(0,a ) - ^(0,a ) - ^^(0,0) = 
Differentiating with respect to y, we get 



(u, + — X — [u, v) 



dydx"^ dy du \ g J ^ dy dv \ g 
The first term of the addition, taken at {Q,u^), is null. For the second term, we 



need to compute the partial derivative of (/4 — /D/ff with respect to v : 

( h-f2 \_K, d . 2N 



dv \ g / g dv ' 

^ h- 3/2/4 + 2/1 /4 - /I % 

o o2 dv 
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Developing the expression of g, we get 

9 = /2/4 - /f /4 - /I - fi + 2/1/2/3 

Let us differentiate with respect to v : 

^ = f2{h-hf2)+f4{h-fi)-fKh-hf2)-2hh{h-hh)-5fi{h-fi) 

- 2/3(/5 - /3/2) + 2/i/2(/5 - /3/2) + 2/2/3(/3 - /1/2) + 2fifs{f4 - /|) 

Taken at (0, 0), each term with even subscript vanishes and we have 

^(0, 0) = cr^(/X6 - ^40-^) + M4(M4 - o-'') - Sa'^ilJ-A - tr^) 
ov 

= 0-^6 - 3/i4Cr^ + 2(7® + (/Z4 - a"*)^ 

Finally 

d_ ( U-fi \ , . ^ M6 - 3(7^/74 + 2^6 _ aVs - 3/^4^^ + 2a« + (^4 - a^)' 

V 5 ' ^ a2(^4-(T4) a4(/X4-(74) 

4 

cr* - /i4 



Therefore 

(0,^') = 0+^(0,^2) - ( ) (0,0) = - — ^ X —J^ = 



dydx^ ' dx ' dv \ g J ' ii^ — a 

This is what we wanted to prove. Let us compute now the fourth partial deri- 
vative of / with respect to x. We have to obtain first an expression of the third 
partial derivative of I with respect to x : 



c^^du^d_(u-f^\, , ^ dv d (h- fi 

dx^ dx du 



, , dv d (fi-fi\ , , 



The only term we do not know is the partial derivative with respect to u of 

(/4 - !l)lg- We have 



ou \ g J g du ^ ' g^ du 



^ h- Uh - 2/2/3 + 2/2^/1 h - fi dg 
9 9^ 9u 

with 

1^ = /2(/5 - /4/1) + /4(/3 - /2/1) - f!{h - /4/1) - 2/4/1 (/2 - /f) 

- 3/|(/3 - /2/1) - 2/3(/4 - /3/1) + 2/l/2(/4 - /3/1) 

+ 2/2/3(/2 - f!) + 2/l/3(/3 - /2/1) 

Notice that this quantity vanishes at (0,0). Therefore the partial derivative of 
(/4 — fi)/g with respect to u, taken at (0, 0), is null as well and we get back 

d^T 
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Differentiating once more, we obtain 



dx^ dx \dx dv? \ g / ' dx dvdu \ g / ' 

+ ^ X — (l±^) (uv) + —X — f A-^) /u v) 

dx"^ du \ g J ^ dx"^ dv \ g / ' 

+ ^ X f ^ X ^ (U,V) + ^ X ^ (l±^) (u,v) 

dx \dx dudv \ g J ' dx dv"^ \ g J ^ 

Let us compute it at (0,cr^) : 



with 



dx^ ' dx \dxdy) ' dx'^dy 

and 

dv? \ g / g dv? ^ g^ dv dv \ g J g"^ dv dv? ^ 

/4-/|92g 2 (dgV 



^(^ g'^ \dv 

Hence 



) (/4-/I) 



\ g ) 

The two remaining terms are the derivatives of quantities which we have aheady 
computed. We evaluate them directly at (0,0), which is straightforward since 
/j (0, 0) = when j is odd : 

|^(/4 - /l)(0, 0) = ^ih- Ufi - 2/2/3 + 2/|/i)(0, 0) = - 3(7^4 + 2a^ 



and 

0(0, 0) - 1^ (1^) (0, 0) = a\,,6 - ^4^2) + M4(M4 - (7^) - - 2m4CT^ 

- 3(7^(/i4 - 0-^) - 2/i| + 2(T^/i4 + 2(7^/14 + 

This is equal to cr^/^e — M4 + 3cr^ — 3/i4cr^ after simplification. Thus we have 
^ (!±lJl\ rn m - ^^Ve - 3aV4 + 2a« - + m1 - + 3m4<7^ 

Finally 



(0,(7 ) = 



j2 



(7° (7° (7° 

We obtain the announced term and the proof is completed. □ 
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8 Around Varadhan's lemma 



We denote by i/n,p the distribution of {Sn/n,Tn/n) under p®" and by 6n,p 
the distribution of {S„/n,Tn/n) under ]ln,p- We saw in section 3 that, if A is 
finite in the neighbourhood of (0,0), then the sequence (^Vi, /))«>! satisfies a 
large deviation principle with speed n, governed by the good rate function I. 
Moreover, for any A cM?, 

. X4nA-exp(^) di^uA^^y) 

Yet we cannot apply Varadhan's lemma directly since A* is not a closed set and 
F : {x, y) I — > / {2y) is not continuous on A. However we have the following 
proposition : 

Proposition 33. Suppose that p is a non- degenerate symmetric probability mea- 
sure on M. such that (0, 0) G Da- Then 

liminf - In / cxp ( ^77- ) di/nA^'V) ^ 

n->+cx) n J ^, \ 2y J 

We assume that there exists r > such that Mr + In p{{0}) < with 

= sup { 1^ : (x, y)GCn BA{(0, 0)} } 

where is the open ball of radius r centered at (0, 0) and C is the closed convex 
hull of { {xjx"^) : X is in the support of p}. If A is a closed subset o/M^ which 
does not contain (0, a^) then 

limsup -In I exp ( ) dunA^^ 2/) < 
n JatiA ^ 2j/ / 

Let us give first some sufiicient conditions to fulfill the hypothesis of the propo- 
sition. 

To ensure that there exists r > sucli that + lnp({0}) < 0, it is enough that 
one of the following conditions is satisfied : 

(a) p has a density 

(b) p({o}) < i/ys 

(c) There exists c > such that p(]0, c[) = 

(d) p is the sum of a finite number of Dirac masses 

Indeed, the function F is bounded by 1/2 on C\{(0, 0)} C A*, thus for any r > 0, 
Mr < 1/2. Therefore, if p has a density, or more generally if p({0}) < e~^/^, 
then for aU r > 0, + lnp({0}) < 0. 

On the other hand, if there exists c > such that ]0,c[ does not intersect the 
support of p (especially if p is the sum of a finite number of Dirac masses) then 

Cc {{x,y) eR^ :c|a;| <y} 
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Therefore 

v(.,rieCnBA((0.0)) g . fli < M < i- 

Hence for any r > 0, < rjlc. Since p is non-degenerate, p({0}) < 1, thus 
there exists r > such that lnp({0}) + r/2c < 0. Therefore the conditions (c) 
and (d) imply that + lnp({0}) < 0. 

Before we prove proposition 33, we need two prehminary lemmas, the following 
one being very useful for handling superior limits. 

Lemma 34. // (ui («.))„>!, . . . , (u/c(^))n>i O'Te. k sequences of non-negative real 
numbers then 

limsup — In ^S^Ui{n) ] = max limsup — In Uj (n) 



We refer to [7] for a proof. The second lemma we need is a variant of the upper 
bound of Varadhan's lemma. Recall that a topological space X is Hausdorff if, 
for any {x, y) € X'^ such that x ^ y, there exist two disjoint neighbourhoods of 
X and y. The Hausdorff space X is regular if, for any closed subset F oi X and 
any x ^ F, there exist two disjoint open subsets Oi and O2 such that F d Oi 
and x e 02- 

Lemma 35. Let X he a regular topological Hausdorff space endowed with its 

Borel a-field B. Let {vn)nei^i c- sequence of probability measures defined on 
{X,B) which satisfies a large deviation principle with speed n, governed by the 
good rate function J. For any bounded continuous function f : X — > R, we 
have for any closed subset A of X, 

limsup - In / e"/^^' d:^„(x') < sup (/(x) - J{x)) 

Proof. Let A, a > 0. We define 

J-\[0,X\) = {x € X : J{x) < X} 

The set J~-'^([0, A]) H A is compact since J is good and A is closed. Since / is 
continuous, J is lower semi-continuous and X is regular, for any x £ A, there 
exists an open neighbourhood 14 of a; such that 

sup f{y) < f{x) + a and inf J{y) > J{x) — a 
yeVa: yeVa, 

The collection Vx, x € J~"'^([0, A]) H A, is an open cover of the compact set 
J~^([0, A]) n A. Let {Vx,, 1 < i < k) he a, finite subcover extracted from this 
covering. Setting U = Uj=i Vx^, we have 

J A Jv,:^ J A\U 

k 

< J2 e"^("')+"«i/„(y,J + e"ll^ll~z.„(^\[/) 
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Moreover 

inf{ J{x) : X € A\U} > inf{ I{x) : x ^ J-^([0, A])} > A 
Therefore, using the large deviation upper bound and lemma 34, 

limsup - In / e"'^'-'"'> dun{y) < max ( max {f{xi) - J{xi) + 2a), ||/||oo - A ) 

n^+oo n J A \l<i<k / 

< max (sup (fix) - J{x)) + 2a, - a) 

We conclude by letting successively a go to and A go to +oo. □ 

Proof of proposition 33. If V is an open neighbourhood of (0,ij^) which is 
included in A* then 

The large deviation principle satisfied by I'n^p implies that 

liminf — In / cxp ( ^77— ) dyn nix, y) > — inf 7 > —7(0, a^) = 

n->+oo n J ^, \ 2y J V 

We prove now the second inequality. Let a > 0. The function I is lower semi- 
continuous on M^, thus there exists an open neighbourhood U of (0,0) such 
that 

y{x, y) eU I{x, y) > 7(0, 0) - a = - In p({0}) - a 

The above equality follows from lemma 26. By hypothesis, there exists r > 
such that Mr + lnp({0}) < thus, by choosing a sufficiently small, we can 
assume that 

Af^ + lnp({0}) + a < 

Since Mr decreases with r, we can take r small enough so that C U. Notice 
next that 

(-,-)=-E(^.,^fc^)eC a.s 
k=i 

therefore, setting C* = C\{(0,0)}, 

/ exp di7„,p(x,y) = / exp ?/) 

Let us decompose 

c* n A c (c* n B^) u (c n B^ n A) 

We have ^ 

/ exp ( ^ ) dvn.p{x, y) < exp(nMr) ?„.p (U) 

The large deviation principle satisfied by Un,p implies that 

1 f (nx'^\ _ 

limsup — In / exp I — — 1 dvn^p{x, y) < Mr — inf 7 < Mr + In p({0}) + a 
„_j.+oo n Jc*nBr V 2?/ / ' k 
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Next, the set CnBJi is closed and does not contain (0, 0) thus the function F 
is continuous on this set. Moreover F is bounded on C*. Hence, by the previous 
lemma, 

limsup-ln / exp (^^^ dun,p{^,y) < sup {F - I) 

n->+cx> n JcnB^nA ^ ^ CnB=nA 

Lemma 34 implies that 

limsup-ln / exp ( ) dun,p{x,y) 
n-)-+oo n JcnA V 2?/ / 

< max ( Mr + In p({0}) + a, sup {F - I) 
\ CnB^nA 

o 

Since p is symmetric and (0,0) G Da, proposition 31 implies that G = I — F 
has a unique minimum at (0, cr^) on A*. Suppose that 

inf G = 

CnB^nA 

Then there exists a sequence {xk,yk)keN in C fl n A c A* such that 
lim G(xk,Vk) = inf G = 

For k large enough, G{xk,yk) < 1/2 thus I{xk,yk) < 1, i-e., {xk,yk) belongs to 
the compact set { {u,v) G : I{u,v) < 1 }. Up to the extraction of a subse- 
quence, wc suppose that {xk,yk)ke^i converges to some (xojj/o), which belongs 
to the closed subset C fl n ^. Moreover G is lower semi-continuous, hence 

= liminf G{xk,yk) > G{xo,yo) > 

Therefore G{xo, yo) = and thus {xq, yo) = (0, cr^) e CflB^flA, which is absurd 
since A does not contain (0,c7^). Thus 



and 



inf G > 

CnBsnA 



max [ Mr + In p({0}) + a, sup {F - I)] <0 

CnB9.nA 



This proves the second inequality. □ 



9 Proof of theorem 1 

The proof of theorem 1 relies on the variant of Varadhan's lemma exposed in 

the previous section. Suppose that p is a symmetric probability measure on K 
with positive variance and such that (0, 0) G Dj^. We assume that one of the 
four conditions given in the paragraph below proposition 33 is satisfied. 
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We denote by 9n^p the distribution of {Sn/n,Tn/n) under /i„,p. We saw in sec- 
tion 3 that for any AcM.'^, 

/a* exp (^) di>n,p{^' y) 
Let U be an open neighbourhood of (0,0-^) in R^. Proposition 33 impUes that 

hmsup — In On^pCU") = hmsup — In / exp I — — J di'n,p{x, y) 

n^+oo n ' n->-+oo U JA'nW^ ^ % ^ 

- Uminf - In / exp ( ) di'n,p{x, y) < 
n^+oo n J ^, \ 2y / 

Hence there exist £ > and no e N such that for any n > no, 



n—^oo 



Thus, for each open neighbourhood U of (0,f7 ), 

lim ]l„^J(^,I^)euA=0 

n-5-+oo \ \ n n / J 

This means that, under /Xn,/9, {Snln,Tn/n) converges in probabihty to (0, cr^). 



10 Proof of theorem 2 

In this section, we first give conditions on p in order to apply theorem 17 to 
the distribution Up. Recall that this theorem states that, for n large enough, 
Vn^p has a density (/„ with respect to the Lebesgue measure on such that, 
for any compact subset Ki of Ai, the admissible domain of /, when n — ^ +00, 
uniformly in {x,y) € Ki, 

5n(x,y)~^(detDf,,^)/)'/%-"^(-.^) 

We use then the Laplace method, as we announced in the heuristics of sec- 
tion 3, to obtain the fluctuations theorem 2. The proof relies on the variant of 
Varadhan's lemma and the expansion of / — in (0, cr^) given in proposition 32. 

o 

We first notice that D\ ^ since it contains Mx] — 00, 0[. Next, we have the 
following lemma : 

Lemma 36. If p has a probability density / with respect to the Lebesgue measure 
on R, then u*^ has the density 

h :ix,y)eR^ -j== f J / J l..<2, 

with respect to the Lebesgue measure on M.^. 
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Proof. Let hhe a bounded continuous function from M? to M. We have 

h{x,y)di^fix,y)= I h{{z,z') + {t,e))dp{z)dp{t) 

= [ h{{z,z^) + {t,e))f{z)f{t)dzdt 

= I h{z + t,z^ +f)f{z)f{t)dzdt+ I h{z + t,z^ + f)f{z)f{t)dzdt 
Jd+ J d- 



I+ + I^ 



with D+ = {{z,t) e M."^ : z > t} a.nd D- = { {z,t) eR^ : z <t}. Indeed, the 
Lebesgue measure of the set { {z,t) G : z = t} is null. 



We define <p : (z.t) gR^ > — > {u,v) = {z + t, z"^ + t"^). If {z,t) G D+ , then 
u'^ _ (z + tV 

T ~ V 2 / 



,2 /._L^\2 z'^+t'^ V 



thus 

ip{D+) C A2 = { {u, v) eR^ <2v} 

Similarly ^{D~) C A2. Conversely, if {u,v) G A2 and u = z + 1, v = z^ + 

then t = u ~ z and 2z'^ — 2uz + {v? — w) = 0. This last quadratic equation in 
z has for discriminant 4(2?; — u^) > 0, thus there are two distinct real-valued 
roots. Therefore 

{z,t) = ]^{u+\/2v-u^,u- \/2v - u^) e £>+ 

or 

[z, t) = -{u- \j2v - V?, u + \j2v- V?) e 

This proves that 9? is a one to one map from £)+ (resp. from D~) onto A2. 
Moreover is on D+ U D~ with Jacobian in [z, t) 



det 



(2. 2O =2k-i|=2v/^^;^^0 



The change of variables given by ip yields 

T T [hi 1 J u + J u- V2^;^ ^ ^ ^ 

'^ = '- = I'^^'^^^Tl^ ^ [ 2 ) f [ 2 ) '^'"^ 

By adding theses two terms, we get the lemma. □ 
Notice that, if p = Af{0, 1), then for any {x,y) G M^, 

This is precisely the formula of proposition 7 for n = 2. 
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By theorem 17, the expansion of (/„ holds as soon as there exists q G [l,+oo[ 
such that /2 e L'^(R'^). However the computation of /2 is not feasible in general. 
Proposition 22 says that the previous condition is satisfied if there exists p € ] 1 , 2] 
such that /2 e L''(]R'') so that the expansion is true. Let us take a look at this : 



\f2{u,v)\^ dudv 

M2 



fP ({u + y/2v-u'^)/2) fP ({u - ^j2v-u^)/2) 

Let us make the change of variables given by 

(u, v) I — > {x,y) = ^ (m + \l2v- u^, u + \/2v- v?) 

which is a C^-diffeomorphism from A2 to (see the proof of the previous 
lemma) with Jacobian in (m, u), 2^/2v — v? = 2{y — x) > Q : 

[ \f2{u,v)\Pdudv= [ 2{y-x)ly:.,dxdy 

(y — X)'' 

By symmetry in x and y, we get 

/ \f2{u,v)\Pdudv= [ fP{x)fP{y)\y-x\^-Pdxdy 

The next proposition follows from theorem 17, proposition 22 and the previous 
equality : 

Proposition 37. Suppose that p has a density f with respect to the Lebesgue 
measure on M such that, for some p £]1,2], 



{x,y)^fPix + y)fP{y)\x 



i-p 



is integrable. Then, for n large enough, i'n,p has a density gn with respect to 
the Lebesgue measure on such that, for any compact subset Kj of Aj, when 
n +00, uniformly over (x, y) G Kj. 



5n(x,y)~^(detD2^,^)/)'^%-"^(-.^) 



Let us prove now theorem 2. Suppose that p is a probability measure on M with 
an even density / such that there exist t;o > and p e ]1, 2] such that 

/ e''°'^f{z)dz< +00 and / fP{x + y)fP{y)\x\^-P dx dy < +00 

The first inequality implies that 

yv<Vo Vm e M ^-voz^ guz+vz^ ^ ^uz-{vo-v)z^ _^ q 

\z\^+cx 

o 

Therefore Kx] — oo,wo[C Da and thus (0,0) G Da- Moreover p is symmetric 
(since / is even) and its support contains at least three points (since p has a 
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density). Proposition 32 implies that / is C°° in a neighbourhood of (0,cr^) and, 
when {x,y) — >■ (0,c7^), 

^ '^^ 2y 2(^*4-^74)^12(78 

where 114 is the fourth moment of p. Denote by B5 the open ball of radius 
6 centered at (0,cr^). It follows from the previous expansion that there exists 
5 > such that for any {x, y) gBs , 

We can reduce S, in order to have c Kj where Kj is a compact subset of 
Ai. Moreover Ai c Dj c A* thus B5 n A* = B5. 

Let n e N and let / : M — > M be a bounded continuous function. We have 

^ An + Bn 



with 



An= I /(a;n^/^)exp d?„,p(a;,?/) 
Bn= I }{xv}''^)exp{^^^ dvn,p{x,y) 

o 

Proposition 37 implies that, for n large enough, v^^p has a density Qn with respect 
to the Lebesgue measure on which verifies, when n — >■ +00, uniformly over 

{x,y) G Ki, 

ff„(x,y)^^(detDf,,^)/)'/%-"^(-'-) 
Let us put e^"^^^'^) in the expression of An : 

= n / /(a;ni/^)e-"°(^'^)i?„(a;,y)dx(i2/ 

where we set il„ = e^^^^''^^ gn(x^y) j n. We define 

B5,„ = { (a;, y) e R2 : a;V\/^ + y V« < } 

Let us make the change of variables given by (a;, y) 1 — )• (a;n~^/*,j/n~^/^ + cr^), 
with Jacobian n~^l^ : 



We check now that we can apply the dominated convergence theorem to this 
integral. The uniform expansion of 5„ means that for any a > 0, there exists 
no e N such that 



(a;, y) G Kj and n>no 



F„(a;,y)27r(detD2^_^)j) - 1 



< a 
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If {x,y) € Bg^n, then {xn,yn) = {xn ^^'^,yn + cr^) e c Kj, thus for all 

n > no and {x,y) € Bg^n, 



Moreover {Xn,yn) — ^ (O)!^^) thus, by continuity, 
Therefore 

\n^'^ yn / n->+c» 
The expansion of G in the neighbourhood of (0, a'^) implies that 

exp {-nG 1,;^. ^ + - J J exp " I^J 

Let us check that the integrand is dominated by an integrable function, which 
is independent of n. The function 

{x,y)^{Bl^^^iy"'^ 

is bounded on by some Ms > 0. The uniform expansion of gn implies that for 
all {x, y) € Bs, Hn{x, y) < Cs for some constant Cs > 0. Finally, the inequality 
(*) above yields that 

iB,,.(.,ri/(x)exp (-„o (-fL,i + ^„ (;^';| + 

and the right term is an integrable function on M^. It follows from the dominated 
convergence theorem that 



By Fubini's theorem, we get 



1/4 / 4\ 

n ' / „, ^ / /i42^ \ 



Let us focus now on i?„. The distribution p is symmetric, it has a density and 
(0, 0) belongs to the interior of Da, thus proposition 33 implies that 

limsup-ln / exp(^^) dvnJx,y) <0 
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Hence there exist £ > and no > 1 such that for any n> no, 

I exp (Iun,p{^,y) < e""^ 

and thus B„ < H/lUe""^ so that B„ = o(nV*). Therefore 

Applying this to / = 1, we get 

where we made the change of variables y = ^^x'^ /{12a^). Finally 

--(/(|7.))r.(St)"'Kl)7/<^'-(-^)- 

An ultimate change of variables s = [i^^xja^ gives us 
This ends the proof of theorem 2. 
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